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The conduction electron density of states nearby single magnetic impurities, as measured re- 
cently by scanning tunneling microscopy (STM), is calculated, taking into account tunneling into 
conduction electron states only. The Kondo effect induces a narrow Fano resonance in the conduc- 
tion electron density of states, while scattering off the d-level generates a weakly energy dependent 
Friedel oscillation. The line shape varies with the distance between STM tip and impurity, in qual- 
itative agreement with experiments, but is very sensitive to details of the band structure. For a 
Co impurity the experimentally observed width and shift of the Kondo resonance are in accordance 
with those obtained from a combination of band structure and strongly correlated calculations. 
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Recently, several groups have demonstrated using 
scanning tunneling microscopy @ | that a magnetic 
Kondo impurity adsorbed on the surface of a normal 
metal causes a narrow, resonance-like structure in the 
electronic surface density of states (DOS), whose asym- 
metric line shape resembles that of a Fano resonance || . 
The experiments were performed with single Ce atoms 
on Ag pf as well as with single Co atoms on Au |^] 
and Cu H surfaces by measuring the I-V characteristics 
of the tunneling current through the tip of a scanning 
tunneling microscope (STM) placed close to the surface 
and at a small distance R from the magnetic atom (see 
Fig. 0(inset)). 

Although the Kondo resonance, formed in the local 
conduction electron density of states (LDOS) at the 
Fermi energy sf due to resonant spin flip scattering, has 
been known for a long time Q, the precise line shape 
was not studied earlier because of the limited spatial res- 
olution available in experiments. Wide tunnel junctions, 
which were proposed as measurement devices Q, probe 
the averaged DOS rather than the LDOS. Such experi- 
ments exhibited a giant zero bias resistance peak |^|8| or 
a weak conductance peak, the latter induced by assisted 
tunneling through impurities in the tunnel barrier . 



In the present Letter we present a detailed microscopic 
study of the Kondo line shape as measured by STM in 
the vicinity of a single magnetic ion. We assume that the 
STM current is predominantly due to tunneling into the 
conduction LDOS, i.e. we neglect direct tunneling into 
the d- or f-level of a Co or Ce ion. This assumption is jus- 
tified, because the d- or f-level is localized deeply in the 
atomic core, and is sufficient to explain the experimental 
findings as seen below. When a discrete (single-particle) 
level is coupled to the conduction electron sea (bare DOS 
po) via a hybridization V, there is a twofold effect: (I) 



the discrete level is broadened and (2) the continuous 
conduction LDOS becomes modified. The resulting line 
shape in the LDOS is called Fano resonance, in reminis- 
cence of the first study of this problem in the context of 
atomic physics [|J. Here we generalize this problem to 
the interacting case, i.e. when the discrete level arises 
from a many-body effect, like the Kondo resonance. It 
is shown that the many-body correlation effects and the 
consecutive Fano line shape in the conduction LDOS may 
be understood in separate steps, thus greatly simplifying 
the theoretical treatment as compared to other studies 
10- 

For concreteness we focus on a Co atom on Au and use 
the Anderson model JlJ] with a fivefold orbital degener- 
acy of the d-level £4 < 0, m = 1 . . . 5, 
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where H Q = e^a^ a a^ a is the kinetic energy of the 
conduction electrons, and a\^. a , d m,a are the electron op- 
erators in the conduction band and in the d-level, re- 
spectively. U is the Coulomb repulsion between two elec- 
trons in any of the local levels. There are two types 
of resonances in the DOS of the impurity d-level: (2) 
the Co d-levels with a broadening A = TrlVppo an d 
(ii) the Kondo resonance [jl3| whose width is given by 
the Kondo temperature Tk- In order to make con- 
tact with experiment, one must obtain realistic esti- 
mates for the parameters of the model. Therefore, we 
have applied the semi-relativistic, screened Korringa- 
Kohn-Rostoker method jl4j in combination with the local 
spin-density approximation (LSDA) ]f5[ | for calculating 
the self-consistent electronic structure of a Co impurity 
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placed onto a Au(lll) surface. The LSDA assumes a 
spontaneous magnetization and consecutive Stoner-like 
splitting of the local level, which is fictitious for a single 
impurity. Although this does not account for properties 
of dynamical origin like the splitting into lower and up- 
per Hubbard states or the Kondo effect, it is an accurate 
method 
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FIG. 1. LSDA result for the LDOS of a Co impurity on 
a Au(lll) surface, showing majority (j) and minority (J.) 
states. The inset shows the experimental setup schematically. 



to determine static quantities like the average impurity 
occupation number n d — J2ma(^m,a-d m a) ■ This is be- 
cause these are determined on time scales much shorter 
than the spin flip time r = h/ksTK, so that the lo- 
cal moment is effectively static for this purpose. Our 
LSDA results are summarized in Fig. [IJ The orbital 
degeneracy is lifted due to crystal field splitting, and 
A ~ 0.2eV. The on-site Coulomb repulsion U is propor- 
tional to the LSDA Stoncr splitting and may be estimated 
as U — 2.8eV fig] . Mainly due to sp-d hybridization, 
the Co d-levels are shifted downward compared to bulk 
Co, so that rid — 8.8 instead of rid — 7 expected from 
the nuclear charge of Co. The excess charge is compen- 
sated by a positive conduction electron depletion cloud 
around the impurity. It follows that of the 11 possible 
charge states of the Co d-level, z = Q, 1, . . . , 10, the sys- 
tem fluctuates only between z = 8, 9, 10. Using Hund's 
rule, these may be identified with the empty (z = 8, 
cr = l, fixed), singly (z = 9, a = ±1/2) and doubly 
(z = 10, a = 0) occupied states of an effective spin- 
1/2 Anderson impurity model without orbital degener- 
acy. The level energy Ed and Coulomb repulsion U of the 
effective model are different from those of the original 
model, Eq. ([!]). They may be extracted as follows Jl7| : 
According to Eq. (Q) and neglecting fluctuations, the 
charge states z = 0, 1, . . . , 10 of the d-level have energies 
E(z) = ze d + Uz(z - l)/2 (Fig. |(a)). Since E(z) must 
have its minimum (ground state) at z = rid — 8.8, we 
have £d/U = —(rid — 1/2). The energies of the singly and 
doubly occupied orbital, measured relative to the empty 



state, E(8), are given by e d = E(9) - E(8), 2e d + U = 
E(10) — E(8) and may be estimated as Ed = -0.84eV, 
U = 2.84eV. The corresponding Co d-orbital spectral 
function A d (u>), calculated using the non-crossing ap- 
proximation (NCA) is shown in Fig. | (b^ In 
addition to the broad peaks located near Ed and £d+U the 
narrow Kondo resonance is clearly seen. The Kondo tem- 
perature, T K = D v / 2A/(ttD) cxp{-l/(2/9 J)}, with_the 
spin exchange coupling J = A/(irpo)[l/\ed\ + l/\Ed + U\], 
may be estimated as Tk — 52K. Here D is a band cut- 
off provided by the bulk and (111) surface [^fj bands of 
Au, respectively; we assume D ~ e f u = 5.5eV. It fol- 
lows from the parabolic form of E(z) (Fig. || (a)) that 
< (>) £d + Ud, if n d < (>) 9. As a consequence of 
level repulsion, the Kondo resonance is somewhat shifted 
upward (present case, n d — 8.8; shown in Fig. || (c)) or 
downward from the Fermi level. 
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FIG. 2. (a) Energy of the relevant charge states z = 8,9, 10 
of the Co ion on a Au(lll) surface, (b), (c) Resulting local 
d-orbital spectral function Ad(ui), T — 4K. The broad peaks 
near Ed and Ed + U and the shift of the Kondo resonance near 
£f are clearly visible. 

According to local Fermi liquid theory |2l] , for T < Tk 
the Kondo resonance is a pure potential scatterer. For 
the analytical treatment below we may, therefore, model 
Ad{oj) — ^Im G d (u — iS) crudely as a sum of three 
Lorentzians. This corresponds to 



G d {uJ - iS) = 



Z d 



lu - Sd-iA uj - £ d -U - 
Zk 

uj — ek — iTj 
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where Ek is the position of the Kondo resonance. Z d , Z\j 
and Zk are the appropriate strengths of the poles. In the 
Kondo regime, at low temperatures (T <g; Tk), unitary 

scattering j2l| implies Z K ~ ttT k /A. 

We now turn to the conduction electron LDOS as mea- 
sured by the STM tip at a distance R from the impu- 
rity. It is related to the electronic field operator smeared 
around the tip position R as "Fr = J ^(x)i7r(x) d 2 x, 
where [/r(x) is a form factor centered at x = R, and 
<F(x) is the two-dimensional Fourier transform of a^ a . 
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According to the Anderson model, the exact conduc- 
tion electron i-matrix is given in terms of the d-electron 
Green's function as t a (iuj n ) — ^■G ( j j(T («a;„). Then 
the correction to the conduction electron Green's func- 
tion due to the presence of the impurity, 5Qfj nCr (iuj n ) — 

GR,a(iUn) ~ SfljrC^n), is 

8&R,a(iUn) = Snjaii^n) ta(iw n ) ^fl^C^n), (3) 

where Q^ a (uJ n ) = -{T T (0)*^))^. The measured 
perturbation in the tunneling LDOS at distance R is 



Using Eq. 



fyflM = -Im 5Q R <T (u - iS). 

TT 

H), it can be expressed as 
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(q R a - l^m^o; - id) + 2q R ^Ret a (uj - iS) 
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where we have defined 
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In the following we drop the spin index a. Im Q^ R (u> — iS) 
and qn depend on R, but on the scale Tk very weakly 
on the energy; thus from here on u) = is taken in these 
quantities. At R = 0, (J_r—o is identical to the asymmetry 
parameter q of the Fano theory Q], and if t(u> — iS) is a 
simple (single-particle) level, t(uj — iS) = (A/ttpo)/(ui — 
Ed — iA), the line shape of Eq. (||) reduces to Fano's well- 
known expression p[oS) = po + dp = po(x + q) 2 / (x 2 + 1), 
with x — (u> — £d)/A [Q. However, Eq. (||) is not limited 
to the non-interacting case. E.g. in the Kondo problem, 
correlations are contained in t(itu n ) and may be treated 
separately (see above), while the Fano line shape is due 
to mixing of real and imaginary parts in Eq. @ . The line 
shape is sensitive both to the scattering phase shift con- 
tained in t(uj — iS) and to the space dependent phase of 
the free conduction electron wave function exhibited by 
q R . In particular, it depends crucially on the local charge 
density and on the details of the band structure . In 

the Kondo regime (\e<i I, [/> A, T<s:T K ), for w~ T K , 
\sk\ < Tk <C A, and using the (approximate) unitarity 

relation Zk ~ ttTk/A (see remark after Eq. (§)) the 
i-matrix may be cast in the form 



Rei(w — i8) = — 



po 1 1 + e l A 1 



Imtfcj — i5) = — I 
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where (3 = -(Z d /ir) e d A/(e 2 d + A 2 ), and e = (w - 
£k)/Tk- In (3 we have taken into account only the 9- 
fold occupied d-state, as this is closest to ef and the 
other charge states give only a small contribution to the 
frequency independent part of t(u>). Thus the final ex- 
pression for the LDOS correction is 



[Im4 0) >-^)] 2 (2q R e + q 2 R -l 
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where Cr — P[2qR — (q R — l)A/ed] arises from potential 
scattering by the d-level and corresponds to a weakly en- 
ergy dependent Friedel oscillation. The first part coming 
from the scattering by the Kondo resonance gives a Fano 
line shape in the tunneling LDOS, controlled by the pa- 
rameter q R . 

Eq. (JsJ) can be fitted to the experimental data for a Co 
atom on a Au (111) surface with excellent agreement 
(Fig. ||). From the fit parameters (see the figure caption 
of Fig. ||) we can conclude that Ek > 0. The value of the 
Kondo temperature Tk — 50K obtained from the fit is 
substantially smaller than the bulk value (Tk > 300 K 
for Co (l3)) as the coupling of the impurity to the con- 
duction electrons is weaker on the surface. Both Tk and 
the shift Ek of the Kondo resonance are consistent with 
the predictions of the NCA calculation in combination 
with the LSDA for Co on a Au(lll) surface. 
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FIG. 3. Fit of Eq. (8) (solid line) to the experimental data 
[2] (squares) at R — 0. The fit parameters are qn=o = 0.66, 
C R=0 = 0.95 and e K = 3.6meV, T K = 5meV~ 50K. 

The remaining task is to calculate q R and C R as a func- 
tion of R. We have considered tunneling of electrons from 
the tip (1) into the 3-dimensional Au bulk states as well 
as (2) into the 2-dimensional Au(lll) surface band p0[ . 
In both cases a free electron-like band structure was as- 
sumed, neglecting, e.g., in case (2) corrections due to the 
herringbone surface reconstruction [M. In case (1) the 



bulk Au Fermi wave number was taken as k 



(b) 



1.21A" 1 , 



while in case (2) it was extracted from the surface band 
structure mapped out in Ref. |2Cj by scanning tunneling 
spectroscopy, fc^ = 0.189A" 1 . A window function was 
adopted for the form factor Ur(x), and the scattering by 
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the Co d-orbital was taken to be isotropic, as the scatter- 
ing phase shift does not strongly depend on angular mo- 
mentum. From Eq. (^) , the Lorentzian line shape with a 
weight A = Re [£^] 2 /( 7r /°o) is formed at those distances 



R where Re^ 0) = {A < 0) or ImC?£ J = (A > 0). For 
case (1), the R dependence of the LDOS is demonstrated 
in Fig. H[ The line shape changes periodically between 
asymmetric Fano and Lorentzian line shape, the period 
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given by the Friedel wave length, A 
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2.6A. 



The overall amplitude decreases with increasing distance. 
For case (2), the results are similar, however with an os- 
cillation period of Xp^ = ir/kp ~ 16. 6A. The latter 
is in good quantitative agreement with the measured ||] 
period of approximately 16A, indicating that the two- 
dimensional Au(lll) surface band gives an important 
contribution to the total LDOS measured by STM. How- 
ever, the precise dependence of the line shape on R is 
not reproduced by our simplifying assumption of a free 
electron band structure. In particular, the model cal- 
culation predicts an initially increasing qn as R grows 
from R = 0, while the fit of Eq. (||) to the experiments 
H implies a decreasing q^. The precise i?-dependence 
of qn will require taking into account the detailed band 
structure as well as the additional scattering phase shift 
induced by the charge of the Co ion and its spatially ex- 
tended screening cloud. This is beyond the scope of the 
present Letter. 




FIG. 4. Qualitative dependence of the line shape of the 
tunneling DOS on the distance of the tip from the impurity 
using the first part of Eq. (^) . 



Summarizing, we have shown that the Kondo reso- 
nance in the magnetic impurity d-level density of states 
(T < Tk) causes a Fano line shape in the density of states 
measured by STM even if we take into account only the 
tunneling into the conduction electron states. As ex- 
pected, the Fano shape arises independently of whether 
the local level is a single-particle orbital or a many-body 
resonance like in the Kondo effect. Superimposed on 
the Fano line is a weakly energy dependent Friedel os- 
cillation induced by potential scattering off the broad d- 
level. From the fit to the experimental data S for Co 
on Au(lll) surfaces we conclude that the Kondo reso- 
nance is shifted upward from the Fermi level, and Tk for 



Co on the surface is substantially reduced compared to 
the bulk value. By combining electronic structure cal- 
culations (LSDA) with methods for strong correlations 
(NCA) these findings are reproduced semi-quantitatively. 
In particular, the shift of the Kondo resonance is due to 
level repulsion between the d-level and the Kondo res- 
onance. We demonstrated the dependence of the line 
shape on the distance of the tip from the impurity by 
using bulk as well as surface state Green's functions. We 
emphasize, however, that details of the band structure 
have to be taken into account in order to reproduce this 
distance dependence quantitatively. 
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